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FURTHER COMMENTS ON THE APPLICATION OF THE METHOD OF 
AVERAGING TO THE STUDY OF THE ROTATIONAL MOTIONS OF A 

TRIAXIAL RIGID BODY, PART 2 


1 . Introduction 

In [A. R., 1971], we described some of the results which we 
have obtained in applying the averaging technique described in 
fF.R.,197l] to the variational equations which arise in treat- 
ing the perturbations of the free rotational motions of a tri- 
axial rigid body. In [A.R.,197l] we carried out the first step 
of the averaging procedure and derived the averaged differential 
equations for a set of canonical variables ( o( , for the 
problem of a triaxial body in a precessing, elliptic orbit 
abput an attracting center. The development was carried out 
to the point that the averaged differential equations [A.R., 
197l^(5»31)] are in a form which can readily be integrated if 
it is so desired. The second step of the averaging procedure 
was not carried out for these canonical variables 
because we were planning to use the averaging technique to 
develop first-order secular solutions for an alternative set 
of noncanonical variables. In [F.R.,1972] we began our dis- 
cussion of the development of these secular solutions by carry- 
ing out the first step of the averaging procedure for a con- 
venient set of noncanonical variables. In the present report, 
we complete the second and final step in the development of 
these first - order secular solutions. 

ft 

References to our earlier reports of June 16, 1970, 

February 19, 1971, August 2, 1971 and February 21, 1972 are 
indicated by [ J.R. , 1971] , fF.R.,197l], [A. R., 1971] and 
fE.R.,1972], respectively. 
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2. First-Order. Secular Solution 

4- j 

Equations (4*24) may be integrated numerically with a 
much longer integration step time than may be used with equations 
(4.3)1 In this Way, we bah obtain the first-order secular solu- 
for the rotational motion under the influence of the 
gravity-gradient torque. We can also, however* integrate the 
averaged system ( 4 * 24 ) analytically with the aid of the inte- 
gral (4.30) . We address ourselves to this problem in the re- 
mainder of this section. 

To begin with, we attempt to integrate (4«24(b)), noting 

9 

that -s^ is .constant. We first express c w (and hence s ,) 
in terms of 6 H bp substituting (b) and ( c) of (4.2£) into 
(4-30) • We find that, to first order. 


2 


V, 


H 


- 1 ( c* 1 h c‘ - a* 1 c 

% 'ET rr 


’ 6 , 


2b 


rr 


0 H + Ch) ’ 


H 


( 2 . 1 ) 


It follows from (2.1) that 


= 


± (A*) l/2 (-c^ 


Q. 


’H 


.#r 3 x- 2 , * -x , 

e H + a 3 °e H + a 2 C e H + A i % + A o ) 


1/2 


, ( 2 . 2 ) 


where 


X 2,2 

A = c> * h 

2 


4b 


1 1 


+ "All equations designated "( 4 . ij) ,i, j, nonnegative integers, 
refer to equations given in [F.R. ,1972]. 
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’ * 4b‘Mi-cj 
A 0 = b_ , 

U 2 2 

c' * h 


A * „ . 8a'*b"C h 
1 2~~Z~ 

C 1 1 h 


* 

a 2 = - 


itb 1 1 2 ( a* 1 

c ,,2 h 2 d* 1 


a c"hC 


« 


h - ij , 


A, = ka . . 

3 i i 

c h. 


If (2.2) is substituted into (4»24(b)), we can write 


dc 




U * 3 * 2 *■ *\ 

-°S H + A 3 c e, +A o°a +a, c Q +A n ) 


1/2 


-(A^^s^dt. (2.3) 

Hr 


+ 


H 


2 e H l % 


In order to integrate (2.3)» we first note that the bi- 
quadratic equation 


k ^3 ^2 

c q + A , c i + A oC 


2"8 h + A l°e H + A o*= 0 


(2.4) 


has the following allowable roots*. 

( i) four distinct real roots* 

(ii) four real roots with two identical roots. 
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(iii) two distinct real roots with a pair of complex roots, 
( iv) two identical real roots and a pair of complex roots. 

In what follows, we consider each Of the four cases. 

Case ( i) s Pour distinct real roots 

In this case, we can write equation (2.3) in the Jacobian 
normal form for real roots [see (BF250.0[|.) and (BF2£0.06) 


* 

g 




-( A ) E 


S). rtt 
0 ° 


(2.5) 


where 



2 , 

172 

f ( a 1 “ a 3^ a 2 -a^) ] 


(2.5) ' 


if a-, , a p , a., and a are the real roots of the biquadratic 

* ^ J T 


equation (4*4-) * and it is assumed that a 1 >a 2 >a^ >a^# 

variable c n is related to u through the equation 
n 


The 


A +A sn^ u^ n * 

Cq = 1 2 , U - ■< U < A. . 

H +A^ sn 2 u* 


( 2 . 6 ) 


t Reference to equations in H and Book of Elliptic Integrals 
for Engineers and Physicists , Byrd,P.F. and Friedman, M.D., 
Springer, Berlin, 1554* are prefixed by the notation BP. 
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Here K is defined by the integral 


(1-k* 2 siri a $ ) 


r' 

where k is the modulus of Jacobian elliptic functions and in- 
tegrals. Explicit values of k^ and A^, i=l,2 ,3>4 will be given 
in the study of the special cases which follows. 

Since the time rate of change of c is real, the value of 

9 H 

c a must either lie between a, and a 0 inclusive or between 
0 H 1 2 3 

and a^ inclusive. We will analyze (2.5) in the subcases which 
follow. 

( 1) a-j>c Q > * The i=l,2>3>4» have the values 

H 


A 1 “ a l^ a 2' a 4 } ' 


A 2 = a^(a 1 - a 2 ) , 


A 3 “ a 2~ a 4 ' 


a 4 = a r a 2 ’ 


(2.7) 


and k A has the value 


k* 2 = <a x' a 2 )(a 3 • V 


( a 1 - a 3 ) ( a 2 - a^) 


( 2 . 8 ) 



6 


for this and the remaining three possiblities under Case ( i) . 
If we integrate (2.5) with respect to u from 0 to u (i.e., 
integrating with respect to time t from t^ to t) > we find 
[see (BF257.00) ] that 


uJ*= + 1_ (A^) l/2 s QQ ^(t-t 1 ), (2.9) 

1 g-y e • L 


where t is the value of t at which c = a n . 

1 9 h 1 

Substituting (2.9) into equation (2.6), we obtain the 

first-order secular solution for c Q and thus the first-order 

n 

secular solution for 6 . 

H 


(2) > c Q > • The i=l,2,3,[|.i have the values 

H 


A 1 

— &2 ( 9-^ — 9 

( a) 

A 1 

= a 3 ^ a l "" a 2 ^ 9 

(b) 

A 3 

* ■* ^ 

( c) 

A 4 

= 8.-^ — • 

(d) 


( 2 . 10 ) 


If we integrate (2.5) with respect to time from - t 2 to t, we 
find [see(BF256.00) ] that 

u? = - 1 (A *) 1 / 2 s o -^- ( t-t 2 ) , (2.11) 

g"* 9 


where t ^ Is the value of t at which 
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Substituting (2.11) into equation 

first-order secular solution for 8^. 

(3) a, < c fl < a 
H 


j. : The A^, i=l,2,3j4* 


(2.6), we obtain the 


have the values 


A 1 - 

a 3 ( % ' 

& 2 ) > 

( a) 

A 2 = 

a 2 (a 3 ' 

y * 

(b) 

> 

u> 

II 

a k - a 2 

9 

( c) 

II 

-3" 

•3 

*3 ' % 

9 

(d) 


( 2 . 12 ) 


If we integrate (2.5) with respect to time t from t^ to t, 
we have [see (BF253«00) ] that 


u 


* 

3 


1 _ 

g* 


(aV /2 




(2.13) 


where t-, is the value of t at which c,, = a-, 

3 e H 3 

Substituting (2.15) into equation (2.6), we obtain the 

first-order secular solution for 8^. 

ri 

(ii) a_ < c n < a )( ; The A* , i=l,2,3>4> have the values 

3 e H 4 


1 = V a i - 

a^) , 

( a) 

2 = a l <a 3 

y ’ 

(b) 

3 = a l " a 3 

9 

( c) 

= a„ - a, 

4 3 4 

9 

(d) 


(2.14) 
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If we integrate (2.5) with respect to time t from t^ to t, 
we find fsee(BF252»00) ] that 

«*= + u*> 1/2 v A (t_ V ’ (2 - ls) 

where t, is the value of t at which c = ai . 

4 % 4 

Substituting (2.15) into equation (2.6), we obtain the 
first-order secular solution for 6 ^. 

The associated first-order secular solution for the var- 
iable "tjJ ^ for each of the subcases described above, is readily 
obtained from (2.1). This eliminates the need to integrate 
( 4 . 24 (a)) directly. Four more integrals of remain to 

be determined. We consider next the variable <b . 

H 

In order to integrate equation (4«24(c)), we first re- 
write it in the form 


x. = (a QO A) _VH - (x,'* - 3x,")cf - x 0 * ' 


s 


0 . 


0 


(2.16) 


H 


H 


2 

We note that if we can integrate c v , /s a and c Q with respect 

r H 2 * * * 6 H 0 H 

to time, to first-order, we can then obtain the solution for <j> • 

H 

From relations (4«25(b)) and (4-30), it is found, to 
first-order, that / 


^H= 


0 


_ h y o 

Z" 


h( c 


1 1 


2 b 


r f 


hfo- 


l' 1 
,r » 


h. 


c h> 


H 


hs 




h -he 


H 


0 


H 


- h ^ 0 


(2.17) 
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If we introduce the identities 


1 = JL ( 1 . + 1 ) , 


h 2 -h 2 2h h+h. 


h-h. 


(a) 


h 2 -h 2 

z° 


1 ( 1 - 1 J > 

^ h-h 2 o h-fh z o 


(b) (2.18) 


a 


.2 v , 2 

h -h z o 


= h ( 1 + 1 )-l , (c) 


2 h+h „ h-h 


equation (2.17) takes the separated form 


’^H = D 1 1 + d 2 1 + 


0 


TT FH 


D. 


H 


’Q. 


h h+c 


H 


0 . 


H 


(2.19) 


where 


D = - 1 (a» 'h 2 - b» 'h + c") , 


D = - 1 (a* *h 2 + b"h + c”) , 


D 3 = a"h . 


We can express the right hand side of (2.19) in terms of 


* 


u if we replace c through the use of (2.6) . We use (2.5) 


to replace t by u* as the variable of integration. We then 


obtain the differential form 
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Vh dt = ? 


9 


H 




r 'M 

D 1 A^-t-A^ sn^u + 


} e° L 


ht A^-A^) 1-^ ^ sn^uT 5 


1 


. _”Z 3 a bn u + E) 

hU 3 +A l’ 1- v 2 sn* u* 3 
o 2 


du 


(2.20) 


where 


-a\ -it 


Ai. - A. 


A 3 “ A 1 


( a) 


( 2 . 20 ) 


2 A, + A- 


(b) 


Upon integrating ( 2 . 20 ) from t^ to t> we find, from (BF336.OI) 
and (BF337.OI) > that 


I li = 


'V H 


dt = - 


’ 9 . 


'H 


+ (if ) i/2 S 0 o X j 


a 3 ( v i (u ?’Ji 2) - v i (0 ’?i 


- \( w !< u r Ti>- V 0 '*?)] 


( 2 . 21 ) 


+ D 2 


h( A3+AX) 


MVvr?- v i (0 ’V?, 


- \{ w i (u f’ tp - w i (0 ’ (Ta’jj - c 3 u i 
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where 1=1* 2, 3 , 4 , and 


v i< u f*r:> 


w i (u i -rj> 


= II(u*y 2 ) (a) 

(j*l,2) (2.22) 


= 1 


|ll(u^, 7^) -F(u* ) 


,(b) 


Here II(u^, y } * j, k positive integer is Legendre* s incom- 
plete integral of the third kind and 7 ^ fa 1 or 7 ^ £ . 

il J 

o o o 

For the special cases where 7 = 1 or 7 = k , the reader 

J J 

may refer to (BFllli06) for appropriate formulas. It can be 
evaluated by using Formulas (BF430) through (BF44 0 ) • 

If, next, we square (2.6), we have that 


c n ~ 

6 h 


ti 

A* ' 1 - y* sn 2 tf' 

k 


( 2 . 23 ) 


where 


- v 2 = a 2 

r 3 a : 


- r = ^ 


(2.24) 


Repeating the procedure used to obtain I , we find that 
(use BF340.02) 
t 


I 2i= 


2 


C/> dt= - 
n 


X- 


S "“" 


a * r 4 

1 1 J V U. + 

3 


+ 2 T3< <4 -y^ K^rrjV 


(2.25) 
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a 


(vr? " V 0 -^ 


( 2 . 2 ?) 


where 


v i ( V r ^ = * 


(a) 


v 2 ( Ul ,y^) = 


2 ( y-- 1, ( ^-^ 




^ E(u i )+(k -f ^ u i 


♦< 2 y“ k**i 2 - “ - * •*• *> 


ii 4 *2, X- 2 

VV 3k 1 


( 2 . 26 ) 


Y j+ sn u A cn u*dri 

^ 2 -V< 

1 ~T i, sn u i 


(to) 


If we take the unperturbed solution ( C / ) H ) 0 as our initial 

value of t the initial value of x will be zero and we can 

H 3 

write 

X rh~ ^H>o“ - (X l” - 3 V )r 2i -V* • 


3 ' 21 3 

(2.27) 

The remaining three integrals of ( 4 . 2 ( 4 .) follow easily. They 
are explicitly 


x = Q* - (Q«) =0 

4 o 

x 5 = <P' - < ^*) 0 =- V’ (t ' 3 I 2 i ) 


(2.28) 


(2.29) 
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(2.30) 


while the unperturbed solutions are taken as initial values 
of the relevant variables. 

Summarizing, we have the secular, first-order solution 
for Case ( i) : 




( c* h 
2b » 



a* 

b* 



+ V 


» 


c = 

e H 


A„ + A an 
1 2 


A^+ A^sh 


~2" 



0 < 



(a) 


(b) 


(2.31) 





(c) 


0' = ( 6' ) o > (d) 

<p ' = ( ^') 0 - x 5 ‘ '(t-3I 21 ) > (e) 

h = (h) q , (f) 

where u^ is given by whichever of ( 2.9) , ( 2.11) ,( 2*13) or 
(2.15) applies to the appropriate subcase. 


Case( ii ) l Pour real roots with identical roots 

This case can be treated as a special case of ( i) and can 
be further grouped into two subcases. 



i4 

{ 1) a-,= a e or a_= a, ^ c_ ; It is seen from (2.8) that 

* 2 e H 3 4 6 h 

k =0. Thus all the elliptic functions reduce to trigonometric 
functions (i.e., sn u = sin u ,etc.) . In either case, equa- 
tion (2.6) becomes 


„ A, +A. 0 s in^u . 
C Q = 1 2 1 , 


H 


A^+A^ s irf.^uj' 


0 < u* < tt /2 , 


( 2 . 32 ) 


where the A^, i=l,2,3>4* are given by either (2.12) or (2.14) 

if a^ = a . They are given by either (2.7) or (2,10) if a = a, , 

Then ^ is given by (2.13) or (2.15) if a ]L = a and it is given 

by (2.9) or (2.11) if as a, , 

3 4 

We note that if we replace (2.3 1(b)) by (2.32) and if, 
for this case (k = 0), we can evaluate 1^ and which 
correspond to (2.21) and (2.25). Then equations (2.31) will 
give us the first-order secular solutions. It is also seen 
from equations (2.21) ,(2.22), (2.25) and (2.26) that if we 
evaluate both II(u*, yj), j=l, 2,3*4 and F(u*) at = 0 then 
these equations will determine bd>th 1^ and I . Formulas 
for the elliptic integrals II{\x^-y 2 ) and F(u^) can be found 
in [1]. Explicitly, they are [see (BF111.01) and (BF121.01) ] . 

F(u*) =u*. 


II ( u i’, rl )ss v lf = 0 


(2.33) 
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tan' 1 ^!- y^) 1//2 tan uj] 

J , if r < 1, 

* j 

(1 "V 


( 2 . 33 ) 


tan 1 [(frl -1) 1//2 tan u^] 2 

J ... 1 * if y . > 1. 

2 172 0 


j 


( 


r j - 1} 


2 


If Y = 1 * we can use (BF111.01) , (BF 121 . 01 ) and (BF111.06), 

J 


and write 


Il(u^il) = tan u^ . 


( 2 . 34 ) 


Thus, the integrals hi and I ^ are determined and equations 
(a) ,(b),(c) ,(d) ,(e) and (f) of (2.31)* together with (2.32), 
give the first-order secular solutions for the six variables 
of interest. 

( 2 ) a 2 = a^i It follows from ( 2 . 8 ) that = 1 . Thus all 

elliptic functions reduce to hyperbolic functions (i.e., 

sn u = tan h u, cn u = sec h u.) . Equation ( 2 . 6 ) takes the form 




+ A2tan h 2 

_________ 

+A, tan h 

4 



0 < 



cO 


( 2 . 35 ) 


where the i=l,2,3,4> are given by either (2.7), (2.10), 
(2.12) or (2.14) and is given by the associated relation 
(2.S), (2.11), (2.13) or (2.15). 
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The analysis proceeds as in the preceding subcase* and 
we have* from (BFlll.O^), that if k* = 1 

= Jln{ tariff - sec <^) * = am Uj* (a) 

II(u£ Y “) 


.[Here y * cannot take on the value one since II(u*,l) = ] . 

Thus the integrals 1^ and I can be determined by (2.21)* 
(2.22) and ( 2.25) *( 2.26) , respectively* and therefore equa- 
tions ( a) * (b) * ( c) , (d) , ( e) * (f ) coupled with (2.35), give the 
first-order secular solutions. 


( 2 . 36 ) 


1 f^n(tan^> i + sec $ 

l -' a i 


(b) 


- v n (i4 


1 +!|/ j V vl/2 

Yj 3 4 


1 *lT'l • (y^D. 


Case ( iii ) l Two distinct real roots and a pair of complex roots 

Let a 1 *a 2 be the real roots and let a , and its complex 

" 5 

conjugate a be the complex roots and assume that a > a . 

1 12 

We can write equation (2.3) in the Jacobian normal form for 
complex roots [see (BP250.05) and (BF250.06) ] 


* A * 
g du = 


dc 


e 


H 


IUl '\ l( V2 )( V3 l( V a 3 l] 


= -(a* : ) 1 ' /2 3 A at. 

m + q° 


(2.37) 
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where g is given by the equation 


* 

g = 


( A* B») 


172 


(23.8) 


and 


A' = f(a^- b ) + a ] , 


B* = f(a 0 - + a “] 


*, 1/2 


a = 


>2 


- 1 (a, - , 

5 3 3 


1 (a + a j . 

2 3 3 


(a) 

(b) 

( c) 

(a) 


(2.39) 


The variable c is now related to u through the equa- 

% 


tion 


a 


'H 


A + A cn u 
= 1 2 , 
Ay- Ayn u* 


0 < u* < 2 k* « 


(2.40) 


^ V. 

where K has the same definition as in Case ( i) and k is the 
new modulus of Jacobian elliptic functions and integrals. 
Since the time rate of change of c is real* the value of 


S H 

c~ must lie between a, and a 0 . If a n > c Q > a 0 * we have 
®H 1 2 1 ®H 2 


A 1 = a l B * + a 2 A * 


( a) 


3 


A. = 


= a^A' - a^B 1 

- A* + B* , 

A* - B* , 


(b) 

( c) 

(d) 


(2.4D 
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and k * 2 has the value 

k** = (a l- a 2 > 2 - (A'-B'l* . ( 2 .^ 2 ) 

1+A* B* 

If we integrate (2.37) with respect to time t from t to t, 
we obtain [see (BF259.00)] 

(A ^) 1 / 2 s 0 A (t-t 2 ) , (2.43) 

v-< 

where t_ is the value of t at which c Q = a 0 . Substituting 

2 e H 2 

( 2 . 43 ) into equation ( 2 . 4 O) we obtain the first-order secular 

solution forcQjj. With the time dependence of c Q known, 

H 

equation ( 2 . 1 ) gives the first-order secular solution for . 

H 

Proceeding as in Case (i), and using equations ( 2 . 40 ) 
and (2.43) in conjunction with equation (2.19), we find that, 
to first-order. 


Y H dts - 

s 

% 




g ..... 
* 072 “ 


>QO 


Si. 


D 1 A 3 


1 +Y cn u 

_3 


h(AyA 1 ) x + ^ cn u * 


D ? A 3 


1 + V cn u 


“1 


h(A^+A^) 1+ y 2 cn u 


+ D 3 


du 


(2.44) 
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where 

y, = \ ~ A g , (a) 

A 3 ' A 1 

y 2 = th 1 A z - (b) ( 24 ?) 

a 3 + a i 

l3 = fk . (c) 

3 

To Integrate equation ( 2 . -t-lj.) with respect to time t , we can 
use (BP361.62) and' rearrange it in the form 


°'^H dt = - ** 

r d a 

1 3 

^ r 3 + 1 " r /Yi \ 


| h(A r A ) | 

L 

^ 1- y^cn u * ) 


D 2 A 3 
h( A^+A^) 


li + 1 ~ r A 2 

r 2 


+ D. 


du 


(2.44) 


1 + Y ^cn u*/ 


If Y ^ ^ 1> and 1# we find, from (BP 36 I. 54 ) [or from 

(BP 341 . 03 ) ], that 


t 




D 1 A 3& + D 2 A 3 Y 3 + D. 
h( Aj-^i Y h( A^+A^) 


u. 
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+ ¥3 




h ( A 3 -V r t 


- R 1 (o,y 1 )J 


+ D A 


y 2 


h("A^+A^5 ^ 2 [ 


R 


x (u ,jr 2 )‘- v o, *V 


, (2.46) 


where 


R. 


!<u , y ± ) = i 


, a -i 

1_ IKu i £ 1 ) - y- f.(u ,-»-*)/ . (a) 

i -y?l rj-i y 


f i( u ^y 2 >=(— 

1 




. 1/2 


tan 


k ^V' a y ?.\l/2 


Li-V 


i -r 


2 / 


sd u 


if y i < k 

Tf-i 


*2 


V - 

sd u , 


if 


*1 

T?-i 


= k 


*2 


(b) (2.47) 


[ /r 1 ) x '\ n 

k* 2 +k*‘ VJ 


1/2 j* 

{.k*a+k*-' a y 2 ^ dn u + 

^2 vi. o 2. l/^ 

( k +k 1 (f ^) dn u - 


1/2 

( T 3 ~li sn u 
172 ■ 

( y -l) sn u 


* 


if ^i > k^ 2 


TT 

V 1 
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2 * , 

If either y^albry^al, the integral in (2*46) is t° 

replaced with the integral given by either (BF^bl.^l) or 

(BF34l*53) i 

Equation (2«4o) can be used to write 
t 




where 



(a) 


(b) 


(2.49) 


Integrating with respect to time* we find from (BF341.O3) and 
(BF341.04) that 




22 


r i<° 


J 4- M' V 2^ 2 -1) -ak » 
\ r 4 / oJ^) ( k**. y *]?•“) 


,*- 2 2.^ ■ 


H.(0, Y > 
4 


+ 2k** R . (u*, 


-i ,n • V H -i <0 ' 


R !< Y h} ~ 


2 X -% 
+ 3n u dn u 
1- y^cn u - * 


- k** r_ 2 (u*. - r_ 2 (0, y , ) 


(2. so) 


where 


TS *7T -1 ^ 

i_l( u > y^) = u + Y4 cos -i (dn u ) 


(2.51) 


R -?( u * y^) S 1 (k* 2 -yj*k ¥ )u*+y, E(u*) (b) 


+ 2 ^ k cos" (dn u ) , 


and y ^ 1 , 

*4 


If y = 1* the integrals in (2 - 14.8) are to be replaced 
with the integrals given by (BF341. 53) and (BF341.54) • Thus 
equations { a ) , (b) , (c) , (d) , (e) , (f) of (2.31) together with (2>4o)» 
(2.46) and (2.50) give the associated first-order secular 


solutions for the six variables. 
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Case ( iv) : Two identical real roots and a pair of complex roots 

This ia a special case of Case (iii) . It can be shown in 

a straightforward manner from ( 2 . 40 ) that if a^ = a^ then 

c = a. and 0 is a constant of the motion. Consequently 
6 H 1 H 

y H> an< * ^ are cons1:ants the motion and <p ^ and 

1 are linear functions of time. 
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